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In this work, we investigate the adsorption of a single cobalt atom (Co) on graphene by means of 
the complete active space self-consistent field approach, additionally corrected by the second-order 
perturbation theory. The local structure of graphene is modeled by a planar hydrocarbon cluster 
(C24H12). Systematic treatment of the electron correlations and the possibility to study excited 
states allow us to reproduce the potential energy curves for different electronic configurations of 
Co. We find that upon approaching the surface, the ground-state configuration of Co undergoes 
several transitions, giving rise to two stable states. The first corresponds to the physisorption of the 
adatom in the high-spin 3d 7 4s 2 (S = 3/2) configuration, while the second results from the chemical 
bonding formed by strong orbital hybridization, leading to the low-spin 3d 9 (S = 1/2) state. Due to 
the instability of the 3d 9 configuration, the adsorption energy of Co is small in both cases and does 
not exceed 0.35 eV. We analyze the obtained results in terms of a simple model Hamiltonian that 
involves Coulomb repulsion (U) and exchange coupling (J) parameters for the 3d shell of Co, which 
we estimate from first-principles calculations. We show that while the exchange interaction remains 
constant upon adsorption (~ 1.1 eV), the Coulomb repulsion significantly reduces for decreasing 
distances (from 5.3 to 2.6±0.2 eV). The screening of U favors higher occupations of the 3d shell 
and thus is largely responsible for the interconfigurational transitions of Co. Finally, we discuss the 
limitations of the approaches that are based on density functional theory with respect to transition 
metal atoms on graphene, and we conclude that a proper account of the electron correlations is 
crucial for the description of adsorption in such systems. 

PACS numbers: 31.15.A-, 31.15.ae, 34.35. +a, 81.05.ue 



I. INTRODUCTION 

Graphene is the focus of many research activities due 
to its remarkable electronic properties* 1 ^ The magnetic 
properties of graphene also attract considerable atten- 
tion because of their potential use in spintronics* 3 - The 
intrinsic magnetism of graphene, being theoretically pre- 
dicted for localized defects or edge-states,— represents a 
challenging problem for experimental studies due to the 
enhanced chemical reactivity of magnetic centers that 
quenches magnetism in the presence of minor contam- 
ination. Furthermore, no magnetic ordering has been 
observed in defective graphene even at liquid helium 
temperatures j 5 -* 6 - which suggests that this kind of mag- 
netism can hardly be employed in practical applications. 
One of the natural ways to induce a magnetic moment 
in a closed-shell compound is to include magnetic impu- 
rities, such as transition metal (TM) atoms, to the sys- 
tem. Apart from being magnetically stable, deposition 
of TM adatoms on the graphene surface may give rise to 
a variety of many-body effects (such as, e.g., the Kondo 
effect)^— known mainly for TM adatoms supported on 
metallic substrates! 10 ' 11 

In contrast to molecular closed-shell adsorbates^ 
light monovalent impurities,— or even some insulat- 
ing substrates ; 14 ' 15 a theoretical investigation of TM 
adatoms on graphene is more challenging due to the pres- 



ence of strong electron correlations. Commonly used 
ab initio concepts, such as density functional theory 
(DFT) or single-determinant Hartree-Fock (HF) -based 
approaches, are often unreliable even in terms of a 
qualitative description of the TM complexes. DFT re- 
sults may vary depending on the parametrization of the 
exchange-correlation functional, whereas the applicabil- 
ity of the HF-based methods is questionable for systems 
with near-degeneracies. The latter is attributed to the in- 
herent limitation of single-reference approaches, such as 
the inability to treat correlation effects arising from the 
mixture of low-lying electronic configurations. Although 
the model Hamiltonian approaches are appealing and 
widely used to investigate strongly correlated systems, 
their results often depend on a number of adjustable pa- 
rameters, which prevent these approaches from being a 
stand-alone tool for studying real materials. 

Among the other first-row transition metals, cobalt 
holds a special place due to the remarkable magnetic 
properties of its compounds related to high Curie tem- 
peratures and large magnetic anisotropics. The experi- 
mentally observed Kondo effect in the system of Co on 
graphene^ followed by considerable efforts in the descrip- 
tion of this phenomenon ; 9 ' 16 ' 17 make Co on graphene one 
of the most attractive candidates for theoretical studies 
among the other TM/graphene systems. 

While there are only a few experimental workS)ii£ 
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there have been numerous theoretical DFT-bascd stud- 
ies that were focused on Co on pristine graphene^i^— 
With respect to chemical bonding, DFT studies predict 
strong covalent interaction arising from the hybridiza- 
tion between the cobalt 3d- and graphenc 7r-orbitals, re- 
sulting in a low-spin adatom configuration (S = 1/2). 
The hybridization leads to relatively large binding en- 
ergies of ~l-2 eV, depending on the initial electronic 
configuration of the Co atom, and also on the particular 
parametrization of the DFT exchange-correlation func- 
tional. Large binding energies obtained within the DFT 
generally imply high migration barriers and, therefore, 
thermodynamic stability of Co atoms on graphene. How- 
ever, this implication seems to be inconsistent with the 
existing experimental evidence of fast surface diffusion of 
transition metal adatoms on graphene, as concluded from 
transmission electron microscopy experiments.— ~— This 
observation suggests a weak bonding, and that there may 
be another mechanism behind the Co-graphene interac- 
tion. 

Not many attempts have been made thus far to go 
beyond the standard local [local density approximation 
(LDA)] or semilocal [generalized gradient approxima- 
tion (GGA)] exchange-correlation approximations with 
respect to Co on graphene. Wehling et al&22- studied this 
system by using the GGA+C/ approach and showed that 
the results depend qualitatively on the on-site Coulomb 
repulsion parameter U. In particular, the application 
of U = 4 eV leads to a rearrangement of the ground- 
state electronic configuration of Co in such a way that 
its spin state changes from S = 1/2 to 3/2 due to the 
3c? — > 4s promotion. Besides, as the U parameter in- 
creases, the adsorption energy of the adatom becomes 
significantly lower, leading to a larger Co-graphene sep- 
aration. Similar results were reported by Chan et al£^ 
The employment of hybrid functionals that incorporate 
some fraction of the exact nonlocal exchange also leads 
to qualitatively similar results, as demonstrated by Jacob 
et al*£- It was also shown in Ref. [l6| that the inclusion 
of the Co-graphene hybridization in terms of a many- 
body theory (one-crossing approximation method) sig- 
nificantly affects the splittings of the 3d cobalt orbitals 
in the crystal field of graphene. All these results indi- 
cate a key role of electronic correlations in the determi- 
nation of the ground-state properties of the Co-graphene 
system, and they suggest that more systematic and re- 
liable approaches should be used in order to properly 
describe electronic properties and chemical bonding of 
Co on graphene. In turn, accurate calculations of these 
properties are necessary to judge how stable Co impuri- 
ties on graphene are, and how the properties of graphene 
itself are affected by the impurities. 

In contrast to previous first-principles studies, in the 
present paper a wave-function-based approach to investi- 
gate a Co adatom on graphene is employed. We use the 
multiconfigurational complete active space self-consistent 
field formalism (CASSCF) along with second-order per- 
turbation theory, which allows us to treat static (strong) 



and dynamical (weak) electron correlations in a balanced 
way. The main advantages this method offers com- 
pared to the DFT are: (i) systematic treatment of the 
exchange-correlation effects, (ii) proper electronic config- 
uration of free TM atoms, (iii) the presence of dispersive 
interactions, (iv) correct asymptotic behavior in the dis- 
sociation limit for open-shell systems, and (v) the pos- 
sibility to study excited states. The price to be paid is 
the fact that the applicability of this method is currently 
limited only to finite and relatively small systems. As 
a result, it is necessary to model the infinite graphene 
lattice by a finite cluster, approximated in this work 
by a coronene molecule (C24H12). This approximation 
looks reasonable if the local adatom-substrate interac- 
tions dominate in the system. Similar cluster approaches 
have been used so far to investigate local correlation ef- 
fects in extended systems by means of the multireference 
quantum chemical methods (see, e.g., Refs. I34T436I) . 



We show that there are two binding mechanisms of Co 
on graphene, which lead to different electron configura- 
tions of Co (high-spin 3d 7 As 2 and low-spin 3d 9 ). The first 
arises from the weak dispersive interaction (physisorp- 
tion) , while the second comes from the strong orbital hy- 
bridization (chemisorption). Despite the chemisorption, 
the adsorption energies are small and do not exceed 0.35 
eV in both cases, which is substantially lower than the 
binding predicted in previously reported studies. We an- 
alyze the results using a simple model Hamiltonian, and 
we conclude that the systematic account of the electron 
correlations is crucial for a proper description of the ad- 
sorption of TM atoms on graphene. Based on this fact, 
we discuss the limitations and applicability of the DFT- 
based approaches to the TM/graphene systems. 



The paper is organized as follows. In Sec. II we present 
the structural model, we briefly describe the CASSCF 
method, and we present other details on the calculations. 
Section III is devoted to the results obtained from first- 
principles calculations, where we first consider a single 
Co atom as a benchmark (Sec. Ill A), and then it de- 
scribes the adsorption of Co on graphene (Sec. Ill B) . The 
nonequivalence of the 3d orbitals of Co and their influence 
on the adsorption are discussed in Sec. Ill C. In Sec. IV, a 
mean-field analysis of the CASSCF results is performed. 
In Sees. IV A and IV B we present a model Hamilto- 
nian and discuss how the parameters of this Hamilto- 
nian can be determined, respectively. In Sec. IV C, a 
discussion on the calculated electron-electron interaction 
parameters is given. The correspondence with previous 
theoretical studies and the limitations of the DFT-based 
approaches with respect to the TM atom adsorption on 
graphene are discussed in Sec. IV D. In the last section 
(Sec. V), we briefly summarize our results and conclude 
this study. 
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II. CALCULATION DETAILS 
A. Structural model 

Many unique physical properties of graphene are deter- 
mined by its gapless cone-shaped band structure, which 
arises as a result of the regular arrangement of carbon 
atoms in the honeycomb lattice. If translational symme- 
try of the lattice is broken, e.g., by contact with defects, 
impurities, or underlying surfaces, the electronic struc- 
ture also changes, affecting the properties of graphene. 
Moreover, graphene should be spatially quite extended 
to retain its properties. In particular, properties of poly- 
cyclic aromatic hydrocarbons (PAHs) are known to con- 
verge extremely slowly with the increase of the number 
of six-membered carbon rings. Theoretical estimations 
show that the highest occupied molecular orbital-lowest 
unoccupied molecular orbital (HOMO-LUMO) gap in 
large PAH molecules closes only at a cluster size of > 
1000 benzene ringsi^l 

However, by investigating the local phenomena, such 
as, for example, interactions with atoms (molecules), an 
exact description of graphene's electronic structure is not 
necessary. Bonding caused by covalent interactions is es- 
sentially local and determined primarily by the atomic 
arrangement as well as by the interatomic interactions 
in the host structure. These interactions vanish very 
rapidly with distance and thus depend only on the local 
surroundings of an impurity. Therefore, it appears rea- 
sonable to consider only a local part of graphene for mod- 
cling this type of interaction. Another aspect of strong 
chemical bonding is the ionic interactions, which scale 
linearly with distance for an unscreened case. In this 
situation, not only do neighboring atoms influence the 
bonding, but a rather large cluster is required to describe 
the interaction in the impurity-substrate system. Ionic 
interactions occur mainly in the case of charge transfer 
between an impurity and a surface. As we will show be- 
low, although in the case of maximal 3d shell filling (3d 9 ) 
some electron transfer from Co to graphene is expected, 
the charge transfer is quite small and any significant role 
of the distant surface atoms can be ruled out. A similar 
reasoning can be applied to the van der Waals (vdW) 
-like interaction due to its weak strength. 

The aforementioned arguments allow us to model the 
surface of graphene by using the finite cluster approach. 
More specifically, we consider a molecule of coronene 
(C24H12) to be a cluster for the local graphene struc- 
ture (Fig. [1} . Similar finite models have been used so far 
in studies on the physisorption of atoms (molecules) on 
graphene (see, e.g., Refs. l38l and l39l). To make the clus- 
ter more realistic, the C-C bond lengths of coronene were 
adjusted to the bond length of graphene (1.42 A). Due to 
symmetry reasons as well as for the sake of computational 
simplicity, we consider only the center of this cluster as 
an adsorption site for the Co adatom. This site corre- 
sponds to a hollow site (77 s coordination) on graphene 
and has been repeatedly predicted to be the most fa- 
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FIG. 1. (Color online) Schematic representation of cobalt 
supported on a coronene molecule used in this work as an 
approximant to the local structure of graphene. 

vorable adsorption site for Co on undistorted graphene, 
according to previous DFT studies. The same adsorption 
site is predicted for the cobalt-coronene complex^ It is 
worth noting that relaxation of graphene may lead to the 
change of the most stable position of cobalt toward a top 
site, as has been shown in the GGA+U studies^^ How- 
ever, the energy difference between hollow and top sites 
was predicted to be quite small. In our study we do not 
take the effects of surface relaxation into account and we 
assume that graphene is locally fiat. 

Due to the Cgv point symmetry, the 3d orbitals of 
cobalt supported on a hollow site of graphene split into a 
single degenerate A\ orbital, and two doubly degenerate 
Ei and Ei orbitals, derived from the d 2 2, d xz {d yz ), and 
d X y (4r 2 -y 2 ) atomic orbitals, respectively. Since the n or- 
bitals of graphene transform according to the E\ and Ei 
representations in the vicinity of the Fermi level, the hy- 
bridization between the cobalt orbitals of A\ symmetry 
{2>d z 2, As) and graphene tt orbitals is suppressed. 

B. Method 

The starting point of an ab initio quantum chemical 
treatment is the HF approximation. The HF method 
is able to describe only the electrostatic and exchange 
contributions to the electron-electron interaction energy, 



whereas correlation effects arc not taken into account, 
-Ejjf = ^totai - -^corr- A hierarchy of widely used post-HF 
approaches [many-body perturbation theories, truncated 
configuration interaction (CI), and coupled cluster (CC) 
methods) 41 - is applicable primarily to the closed-shell sys- 
tems, where correlation effects are rather weak and do 
not significantly affect the electronic structure of the sys- 
tem under investigation. This type of correlations is usu- 
ally referred to as dynamical correlations* 4 ^ according to 
the nomenclature commonly used in the quantum chem- 
ical literature. On the other hand, strong correlations 
may substantially affect electronic structure, so that the 
initial HF wave function (single Slater determinant) be- 
comes unreliable and not suitable for further corrections. 
This is often the case for open-shell systems with near- 
degeneracies, including TM compounds. Effects that 
cannot be properly described at the single-determinant 
level are usually referred to as static correlations. Total 
correlation energy can thus be given as the sum of the 
static and dynamical terms, -E C orr = E sta t + -Edyn' 43144 

One of the efficient ab initio approaches that allows 
systematic treatment of the static correlations is the 
CASSCF method^ In this method, the wave function 
is constructed full CI expansion within a limited 

set of active orbitals, for which all possible occupations 
are allowed. The other orbitals (inactive) have occupa- 
tion numbers exactly equal to 2 and are treated at the 
HF level. Importantly, not only expansion coefficients 
but also the orbitals are simultaneously optimized in the 
CASSCF, which makes this method quite flexible. Dy- 
namical correlations can be taken into account on top 
of the CASSCF wave function using multireference per- 
turbation theories. In this work, we employ the nonrela- 
tivistic CASSCF method in conjunction with the strongly 
contracted variant of the second-order TV-electron valence 
state perturbation theory (NEVPT2) 4 £ as implemented 
in the ORCA program packaged 

In all calculations we employed the Ahlrichs triple- C 
valence basis set (TZV)^ for all atoms. The cobalt atom 
and its six neighboring carbon atoms (central ring in 
the coronene molecule) were additionally supplemented 
by (2df)- and (d)-polarization functions (TZVPP and 
TZVP basis sets), respectively^ We did not apply any 
basis set superposition error (BSSE) corrections to the 
calculated energy values. The counterpoise correction 
(CP)j££ widely used in highly accurate quantum chem- 
ical calculations, does not usually lead to significant 
changes in the interaction energies for weakly bounded 
systems, whereas its application to strongly bounded 
systems is controversial^^ Moreover, the reliability of 
the CP with respect to the multiconfiguration meth- 
ods appears questionable. As a starting point for the 
CASSCF calculations, we used the natural orbitals ob- 
tained at the level of unrestricted single-reference M0ller- 
Plesset second-order perturbation theory (MP2). The 
calculated energies and orbital gradients were converged 
within 10~ 8 and 10 -4 a.u., respectively. It should also be 
noted that upon evaluation of four-index integrals, we did 
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FIG. 2. (Color online) The set of orbitals included in the 
active space within the CASSCF approach (in the form of 
natural orbitals). The symmetry of the orbitals is shown in 
parentheses. 



not use any simplifications such as, e.g., the resolution of 
the identity approximation (RI)j^ 



C. Active space 

The active space used in this work for the CASSCF 
calculations consists of five 3d Co orbitals, one 4s Co 
orbital, and two pairs of bonding (7Tf,) and antibonding 
{n* b ) orbitals of the coronene molecule. This gives rise to 
13 correlated electrons in 10 active orbitals [CAS(13,10)]. 
In Fig. [2l we show schematically the active space in the 
form of natural orbitals obtained at the CAS(13,10) level 
for the ground state of Co at the distance of d = 3.1 A 
from the surface. 

By selecting the surface orbitals, we were guided by 
the following criteria: (i) the orbitals should be localized 
in the center of the surface, i.e., in spatial proximity to 
the Co atom; (ii) the orbitals should be either of E\ or 
Ei symmetry, since the contribution of the orbitals of A\ 
symmetry is not likely in the graphene-impurity bonding; 
(iii) the orbitals should be close to the Fermi level and 
represent a bonding- antibonding pair in order to mimic 
the 7r and 7r* orbitals of graphene. The chosen surface 
orbitals obey these criteria. Moreover, except for the A\ 
metal orbital, for all the other metal orbitals there is 
a corresponding ligand pair of the same symmetry. It is 
worth mentioning that the shape of the selected 7r orbitals 
of coronene is similar to the maximally localized Wannier 
functions associated with the 7r orbitals of graphene^ 

We note that in order to avoid unreasonably expensive 
calculations, our choice of active orbitals does not include 
explicitly the 3d double shell effect, which is still needed 
for an accurate treatment of the excitation energies in 
TM atoms' As we will show, however, the size of the 
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chosen active space is sufficient to obtain a reasonable 
agreement between calculated and experimental values 
of the 3d -O- 4s promotion energies for a single Co atom. 

III. FIRST-PRINCIPLES RESULTS 
A. Single cobalt atom 

The filling of the d-shells of TM atoms is known to 
occur in competition between the d and s states. As a 
result, low-lying excitations in the spectrum of a single 
atom are composed exclusively of the 3d ■(-> 4s transi- 
tions. This competition is more pronounced upon forma- 
tion of TM complexes where TM-ligand bonding takes 
place. While all first-row TM atoms accommodate elec- 
trons in the 4s shell, the occupation of this shell becomes 
much less favorable upon contact with a ligand due to the 
Pauli repulsion between the ligand and the very diffuse 4s 
orbital of the TM atom. Therefore, for investigating the 
adsorption of TM atoms, it is of great importance to take 
into account the interplay between different electronic 
states. Besides the general consideration, the important 
role of the 3d and 4s orbitals in the interaction of Co 
with graphene has been shown by previous studies.— i 16 i 28 

To investigate the role of different electronic states in 
adsorption of Co on graphene, we consider five differ- 
ent electronic configurations of Co: (a) 3d 7 4s 2 (5=3/2), 
(b) 3d 8 4s! (5=3/2), (c) 3d 7 4s 2 (5=1/2), (d) M 8 ^ 1 
(5=1/2), and (c) 3d 9 4s° (5=1/2). As one can sec, these 
states correspond to two quartet and three doublet states 
with different occupancies of the 3d and 4s Co shells. The 
experimental ground state of a single cobalt atom corre- 
sponds to a high-spin quartet state, 3d 7 4s 2 . Unlike the 
LDA/GGA results obtained within the DFT^ this state 
can be predicted even at the Hartree-Fock level of the- 
ory. However, in order to obtain reasonable excitation 
energies, a multireference treatment is required.— 

In Table U we show the excitation energies of a single 
Co atom calculated at the NEVPT2 level for the consid- 
ered electronic states. These states can be assigned to 
the lowest atomic terms of the corresponding configura- 
tions in the experimental spectrum as follows (in ascend- 
ing order): a A F (ground), b A F , a 2 F, a 2 G, and c 2 D^ 
The energies of these terms, averaged over the spin-orbit 
components mj, are given in Table|U The states b 4 F and 
a 2 F correspond to the two low-lying excited states of Co, 
while the others (a 2 G and c 2 D) relate to much higher ex- 
citations. One can see that the calculated energies of the 
first two excited states are in very good agreement with 
experiment, while for the third excited state the agree- 
ment is slightly worse. The maximum discrepancy can 
be seen in the case of the highest energy level. This is 
not surprising in view of the fact that the active space 
used in this work is restricted only to the 3d and 4s or- 
bitals, whereas higher orbitals arc necessary to properly 
describe the high-energy excitations. Another necessary 
condition for an accurate prediction of the excited cn- 



TABLE I. Calculated and experimental excitation energies of 
a single Co atom corresponding to transitions from the ground 
(3d 7 4s 2 , 5=3/2) to low-lying excited states associated with 
the electronic configurations considered in this work. Exper- 
imental energies correspond to the m.j averaged value of the 
lowest atomic term of a given configuration (Ref. lo8h. 
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ergy levels is rclativistic corrections^ which were also 
not utilized in the present work. Nevertheless, despite 
some quantitative discrepancies, the agreement obtained 
between theoretical and experimental spectra for a single 
Co atom can be considered as satisfactory. The consis- 
tency achieved with experimental data shows the reliabil- 
ity of the computational approach as well as an adequate 
choice of parameters for studying the Co-graphene sys- 
tem. 



B. Adsorption of cobalt 

To reproduce the potential energy curves, we perform 
a series of state-specific calculations for different separa- 
tions between Co and graphene. We define the adsorp- 
tion energy of Co in the state i as 

E lds( d ) = E lmf+Co( d ) ~ E 8urf ~ ^Co° Und ' i 1 ) 

where El ur{+Co (d) * s the energy of the interacting system 
in the i state of Co separated by the distance d from 
the surface, -E sur f is the energy of the clean surface, and 
E^° und * s * ne energy of the isolated Co atom in its ground 
state. Equilibrium adsorption energies for different Co 
states can be obtained by taking the minimum of the 
corresponding adsorption energy functions, £ads(^eq) = 
min[£'* ds (d)]. Since the NEVPT2 approach employed in 
this work is size-consistent, the dissociation limit can be 
correctly reproduced, and at large distances Eq. (JXJ) gives 
the interconfigurational energies between the ground and 
excited states of a single Co atom, El ds (d — > oo) = E % Co — 

^ground 
^Co _ • 

In Fig. [3j we show the lowest adsorption energy curves 
calculated at the NEVPT2 level for different configura- 
tions of the Co adatom supported on graphene. The 
adsorption energies, equilibrium distances, and orbital 
occupations for each configuration are summarized in Ta- 
ble HIl As expected, the shapes of the curves as well as 
equilibrium distances do not depend significantly on the 
spin and depend mainly on the orbital occupancies. In 
the case of the same occupations, the low-spin curves lie 
higher in energy due to the difference in exchange energy. 
In the short-distance range, this difference gradually de- 
creases, which can be associated with the loss of the d — d 
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FIG. 3. (Color online) Potential energy curves for different 
electronic configurations of Co on graphene. Zero energy cor- 
responds to the ground-state energy of the noninteracting sys- 
tem. 



exchange energy effect that is known for molecular sys- 
tems containing TM atoms^ 

The adsorption of transition metals differs from the 
adsorption of monovalent or closed-shell adsorbates. As 
can be seen from Fig. [31 the ground-state electronic struc- 
ture of Co undergoes several transitions upon adsorption. 
In turn, the adsorption energy curve (the minimum en- 
ergy curve) is composed of three different potential en- 
ergy curves. As the atom approaches the surface, its 
most favorable configuration (quartet, 3c? 7 4s 2 ) changes to 
3d 8 4s 1 (quartet) at a distance of ^2.6 A, and to 3g? 9 4s° 
(doublet) at shorter distances (~1.8 A). This behavior 
can be qualitatively understood by noting that in the 
vicinity of the ligand, the occupation of the 4s orbital 
becomes less favorable due to the Pauli repulsion, while 
the occupation of the 3d orbital becomes more favorable 
due to the screening of the Coulomb repulsion in the 3d 
shell of Co. We analyze these effects in more detail in 
Sec. IV C. The energies of the two other doublet con- 
figurations lie higher over the entire range of distances. 
The absolute values of the adsorption energies are quite 
small, reaching 0.32 eV in the high-spin 3d 7 4s 2 case, 
while in the other cases the values are even smaller (see 
Table |TT| . These values are significantly smaller than the 
adsorption energy predicted by the DFT studies (~l-2 
(AO^— — — — — — but are consistent with experimen- 
tal observations of the fast surface diffusion of metal 
atoms on graphene^ which suggest weak bonding in the 
system. 

Although the minimum energy curve shown in Fig. [3] 
exhibits three extrema, only two of them (low-spin, 
3d 9 4s° and high-spin, 3g? 7 4s 2 ) are well separated and can 
be regarded as stable configurations. The third mini- 
mum (high-spin, 3d 8 4s 1 ) corresponds to the metastable 
state. Among the two stable configurations, the high- 



spin minimum lies lower in energy than the low-spin one 
and, therefore, represents the global minimum. As the 
low-spin configuration is less favorable in absolute val- 
ues, there is an energy barrier to reach this state. On 
the other hand, if the adatom is trapped in the low- 
spin state, another barrier has to be overcome to re- 
turn in the high-spin configuration. The energy barriers 
can be estimated as a difference between the transition 
state and high-spin (low-spin) state energies. As can be 
seen from Fig. [3J there is a transition state at a distance 
of 1.85 A with an energy of 0.24 eV. Taking this state 
into account, we obtain the energy barriers of 0.56 and 
0.29 eV, respectively for the high-spin— ^low-spin and low- 
spin— ^high-spin transitions. We note that a qualitatively 
similar multiple- minima structure of the potential energy 
curves has been reported for graphene adsorbed on metal 
surfaces.— ~— 

Let us analyze the nature of the Co-graphene interac- 
tions. The maximum adsorption energy obtained for Co 
(0.32 eV) is comparable with adsorption energies typical 
for small closed-shell molecules adsorbed on graphene, 12 
where only weak vdW interaction plays a dominant role. 
In fact, the same mechanism is responsible for the ad- 
sorption of Co in the 3d 7 4s 2 and 3d 8 4s 1 configurations. 
In these cases, the hybridization between the localized 3c? 
orbital of Co and the it orbital of graphene is suppressed 
due to the Pauli repulsion between graphene and the dif- 
fuse 4s Co orbital. Interestingly, however, the surface 
repulsion favors some hybridization between the 3d,A 1 
and 4s Co orbitals in the low-spin 3g! 8 4s 1 state, as in- 
dicated in Table HT1 In the cases where the 4s orbital is 
occupied, the absence of the hybridization with the sur- 
face can be concluded by analyzing the occupancies of 
orbitals associated with the 3c? Co orbital. Lowdin pop- 
ulation analysis^ as well as the occupation numbers of 
the natural orbitals^ show that at equilibrium distances, 
the Co orbitals do not have contributions from/to the 7r 
orbitals of graphene, so that the corresponding occupa- 
tion numbers are close to integer values (see Table [TTJ) . 
It is also worth noting that the binding of Co with the 
nonempty 4s orbital occurs only at the level of second- 
order perturbation theory and not in the pure CASSCF 
calculations, which also confirms that the interaction be- 
tween Co with a nonempty 4s shell and graphene is of 
the weak vdW-type. 

The situation with the low-spin configuration (3d 9 4s°) 
is different. In this case, as can be seen in Fig. [3l the 
potential energy curve exhibits a marked linear depen- 
dence within the distance range of 2.8-4.7 A. This be- 
havior suggests an ioniclikc interaction between graphene 
and Co in the 3g? 9 4s° state. Lowdin population analy- 
sis shows that in this range, the occupation of orbitals 
changes from 3c? 8 85 to 3c? 8 ", which means that some 
fraction of Co electrons is transferred toward graphene. 
This can be understood in terms of the ionization en- 
ergy lowering associated with the filling of the 3c? orbital. 
Therefore, the observed ioniclike interaction can be ex- 
plained as a result of a minor charge doping. At shorter 
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TABLE II. Calculated adsorption energies, equilibrium distances, and orbital occupations for Co on graphene in different 
electronic configurations. Occupation numbers correspond to the equilibrium distances and are given in terms of Lowdin's 
charges^ and natural orbital occupations.— 





S= 


=3/2 




5=1/2 






3d' 4s^ 


3d 8 4s 1 


3d'4s^ 




3d a 4s u 


-Bads, eV 


-0.32 


-0.17 


1.89 


0.24 


-0.05 


d c q, A 


3.1 


2.5 


3.1 


2.3 


1.6 


Orbital occupations^ 


3d 7.00 4s 1.89 


3d 7.98 4s 0.87 


3d 7.00 4s 1.89 


3d 7.98 4s 0.86 


3d 8.22 4s 0.06 


Orbital occupations_ 


3d 7 4s 20 


3d 8 4s 10 


3d 70 4s 20 




3d 8 - 7 4s 00 



a Based on Lowdin's population analysis— 
b Based on natural orbital analysis— 

c In this case, the 3d J \ l Co orbital slightly hybridizes with the 4s orbital forming a 3d— 4s hybrid, (sd)± = l/v / 2(3dA 1 i 4s). 



distances, the 3d 9 4s° adsorption curve exhibits a deep 
minimum typical for strong covalent interactions. In this 
situation, the hybridization between the 3c? Co and it sur- 
face orbitals is more likely since the 4s Co orbital is empty 
and the atom-surface repulsion is much weaker. As can 
be seen from Tabic [Til the natural orbitals of Co in the 
3d 9 4s° state have noninteger occupation numbers, which 
suggests that there is an additional, many-body contribu- 
tion to the Co-graphene bonding. Indeed, analysis of the 
wave function shows that the weight of the principal con- 
tribution in the CASSCF expansion, \^b^A 1 ^E 1 < ^E 2 7r 2b) ^ 
is only 0.84. The remaining part of the wave function is 
distributed over many configurations, which mainly cor- 
responds to the hopping of electrons from the 3d# 2 Co or- 
bital to the ir* b surface orbital. Only one of these config- 
urations reaches a contribution larger than 1%, namely, 
Hd^dl^d^lKl) (4%)- This allows us to con- 
elude that the deep minimum observed in Fig. [3] for the 
adsorption of Co in the 3d 9 4s° configuration is attributed 
to the donation of the metal electrons to the ligand. The 
donation effect can be more clearly seen in Fig. @J where 
we show isosurfaces of the natural orbitals associated 
with the Co 3dg 2 and ir* b surface orbitals. Since the 
donation process involves orbitals of specific symmetry 
(3<1e 2 ); the degeneracy of the 3d Co orbitals is broken 
and an energy splitting is expected between different 3d 
orbitals. We discuss this issue in the next section. It 
is also important to emphasize that the effect of con- 
siderable energy lowering of the 3d 9 4s° Co state has an 
essential multiconfigurational character and thus cannot 
be properly described at the single-determinant level. 



C. 3d-orbital splitting 

Up to now, we have not paid any attention to the 
splitting of the 3d-orbitals, assuming that the electrons 
occupy the lowest energy states in a given configura- 
tion. It is clear, however, that when a TM atom ap- 
proaches the surface, the crystal field tends to split the 
3d-orbitals more strongly the smaller the distance. More- 
over, symmetry-specific hybridization discussed in the 
previous section suggests that these effects are not neg- 




FIG. 4. (Color online) Isosurfaces of the natural orbitals as- 
sociated with the hybridized de 2 (adatom) and 7r„ b (surface) 
orbitals. The numbers correspond to the natural orbital oc- 
cupation numbers. 

ligible. As we already mentioned, if a TM atom is de- 
posited on a hollow site of graphene, its otherwise degen- 
erate 3d-orbitals split into three groups according to the 
A\, Ei, and E2 irreducible representations of the Cq v 
point symmetry. In what follows, we assume that the 
3d-orbital splitting does not depend on the particular Co 
state, i.e., the same for different electronic configurations. 

To estimate the 3d-orbital splitting, we perform ad- 
ditional energy calculations for the 3d 9 4s° configura- 
tion of Co, which correspond to different occupations of 
symmetry nonequivalent orbitals (dx,, d^, and d_E 2 ). 
The potential energy curves for these states are given in 
Fig.G^a). One can see that at distances larger than 2.5 A, 
the energies are almost identical, i.e., they do not depend 
on the particular occupation of the 3d orbital. As the 
separation between Co and the surface decreases, differ- 
ent orbital occupations result in different energies. This 
reflects the fact that the Co orbitals interact differently 
with the surface, which is also clear from the hybridiza- 
tion between the 3ds 2 Co and ir* b surface orbitals, as 
discussed above. As a result of hybridization, the states 
with fully occupied 3d£ 2 orbitals are much more favor- 
able, as can be seen from Fig. [5ja). 

Energy differences between the states of different 3d 
orbital occupations can be interpreted in terms of the 
simple crystal-field picture. If we express the energy of a 
particular state as a sum of one-particle energies, 

£ = ^e^, (2) 
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the effective crystai-ficld splitting parameters, A !3 - = 
Si — ij, can be easily estimated from the total energies 
of different orbital occupations. In Fig. EJb), we show 
these parameters as functions of distance. First, it should 
be noted that the orbital energy levels follow the order 
£e 2 < < £_Ei, which does not depend on the dis- 
tance. This order is consistent with previous DFT— £2- 
and many-body^ studies. As can be seen from Fig.[SJb), 
the Aa 1 -e 1 splitting parameter is relatively small, not 
exceeding 0.5 eV even at short distances. On the other 
hand, the splitting between the E 2 and E\ (Ai) orbitals 
is significantly larger, reaching 4 cV at a distance of 1.6 

A. 

There are two main reasons that can account for the 
large splitting. The first arises from the direct interac- 
tion of Co electrons with the surface, resulting in the hy- 
bridization and/or interorbital repulsion. The second can 
be attributed to the screening effects associated mainly 
with the reduction of the Coulomb repulsion in the Co 
3d-shell. Indeed, as we have shown before, the 3dE 2 Co 
orbital hybridizes with the surface, lowering the energy 
of this orbital. At the same time the hybridization causes 
derealization of the electrons occupying this orbital (see 
Fig. 01), which in turn reduces the Coulomb repulsion 
and leads to further energy lowering. In the next sec- 
tion, we analyze these effects quantitatively based on a 
simple model Hamiltonian. 

IV. MEAN-FIELD ANALYSIS 

A. Model Hamiltonian 

Analysis of a many-body wave function in terms of the 
usual single-particle quantities is challenging. Although 
the single-particle representation of many-body charac- 
teristics can in principle be constructed (e.g., in the 
many-body Green's function formalism), the CASSCF 
wave function cannot be readily mapped onto the single- 
particle picture. As a result, a number of physically 
relevant quantities, such as, for example, the density of 
states, are not available in the multirefcrence methods. 
In this work, in order to more deeply understand the 
physics behind the adsorption of Co on graphene, we 
map the CASSCF results onto a mean-field model. To 
this end, we consider a simple Hamiltonian: 

i<7 2,(7 

Here i and a label orbital indices and spin-projections 
for 3d and 4s electrons of Co (a implies spin inver- 
sion), respectively, £j is the energy of the i electron in 
the absence of electron-electron interaction, and is 
the particle number operator. The parameters Ui, U^, 




1.6 1.8 2 2.2 2.4 2.6 

Distance (A) 

FIG. 5. (Color online) (a) Potential energy curves for different 
orbital occupations of Co in the 3d 9 configuration; (b) energy 
differences between nonequivalent 3d-orbital levels of Co (the 
effective crystal-field splitting parameters) calculated as func- 
tions of the adatom-surface separation (see text for details). 
The splitting of the energy levels is shown schematically inside 
the lower panel. 



and Jij determine intraorbital and interorbital Coulomb 
interactions as well as interorbital exchange interaction 
(Hund's rule coupling), respectively. The first term in 
this Hamiltonian describes the one-electron contribution 
to the energy, including the interaction with the crystal 
field of the ligand. The rest represents the Coulomb inter- 
action Hamiltonian in the Kanamori parametrization^ 
which describes the electron-electron interactions. In this 
paramctrization, the Coulomb and exchange parameters 
for the 3d electrons satisfy the relationship U — U' = 2 J 
if the system is rotationally invariant. It should be noted 
that although the given Hamiltonian does not explicitly 
contain atom-ligand interactions, these interactions are 
implicitly taken into account in the first term through the 
crystal-field splittings of the Co orbitals. In turn, since 
the U, U', and J parameters are defined only for the Co 
part of the system, they should be considered as effec- 
tive. For the sake of simplicity, we also assume that the 
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rotational invariance condition holds and the electron- 
electron interaction parameters have the same values for 
different 3c? orbitals (i.e., averaged over orbitals). Fi- 
nally, we neglect the intraorbital Coulomb repulsion be- 
tween electrons in the 4s orbital (U^s), and interorbital 
3d-4s repulsions (U' s _ d ), taking into account the delocal- 
ized character of the 4s orbital. The interaction between 
the 4s and 3d electrons is taken into account through the 
exchange interaction parameter J s -d- 



B. Determination of model parameters 

Since the Hamiltonian given by Eq. ([3]) is diagonal, the 
energies of different atomic states of Co can be readily 
expressed in terms of the parameters Ud-d, U' d _ d , Jd-d, 
and Jd-s, and the 3g?-4s orbital energies £Ai, EEd £,e 2 j 
and £4 S . By comparing the energies of different states 
obtained in the model with the values obtained from 
the first-principles calculations, the magnitude of the 
electron-electron interaction parameters can be extracted 



-J! 



if the crystal-field splitting parameters, Ay = e 
are known (see the Appendix for the specific expres- 
sions). Unfortunately, these parameters do not gener- 
ally correspond to the effective Ay parameters obtained 
in Sec. HI C. As we already discussed, apart from the 
direct interactions with the crystal field, the Ay param- 
eters may contain screening effects, enhancing the split- 
ting due to additional lowering of the 3g?£ 2 energy level. 
However, there is no straightforward way to establish a 
precise quantitative relationship between these two ef- 
fects from first-principles calculations. 

The splitting in metal atoms due to the presence of 
neutral ligands is in general not large, especially for low- 
coordinated complexes such as atoms on surfaces. Taking 
this fact into account, it appears reasonable to employ 
significantly smaller splittings than Ay upon determina- 
tion of the parameters of the model Hamiltonian. To be 
more specific, we scale the effective splitting by a factor of 
0.2, so that Ay = 0.2 Ay. At the distance of d = 1.6 A, 
this choice is close to the splitting obtained by means of 
the standard DFT-GGA calculations without using the 
[/-corrections (Ae 2 -E! = 0.8 eV, A^-^ = 0.5 eV)£ 
The use of smaller scaling factors does not significantly 
affect the results, while the larger factors lead to an un- 
physical (nonmonotonic) dependence of Ud-d with dis- 
tance. This can be considered as yet another criterion 
for the physically reasonable choice of the scaling factor. 
Finally, we note that the uncertainty in the determina- 
tion of the crystal-field splitting affects only the Coulomb 
repulsion parameters in the short-distance range. 



C. Analysis of the electron-electron interactions 

In Fig. |51 we show the calculated electron-electron in- 
teraction parameters as functions of the Co-graphenc dis- 
tance. As one can see, the exchange coupling parameters 




1.6 1. 



2.2 2.4 2.6 2.: 
Distance (A) 



FIG. 6. (Color online) Electron-electron interaction parame- 
ters for Co on graphene as a function of distance: exchange 
couplings (J s -d and Jd-d) and intraorbital Coulomb repulsion 
(Ud-d)- The latter is shown for different splitting parameters 
Ay (see text for details). S^J^ quantifies the Pauli repulsion 
between the 4s (Co) and ivt (surface) orbitals.— The dashed 
lines mark the most favorable electronic configuration of Co 
for a given range of distances (in accordance with potential 
energy curves given by Fig. [3J). 



(J s -d and Jd-d) are almost constant and do not change 
over the entire range of distances. Jd-d has a typical 
value commonly used in the LDA+J7, and model Hamil- 
tonian calculations of TM compounds^ and it amounts 
to ~1.1 eV. As expected, the s-d exchange interaction is 
significantly smaller (0.1-0.2 eV) and does not play an 
important role in determining the electronic structure. 
In contrast, the intraorbital Coulomb repulsion (Ud-d) 
significantly decreases for shorter distances (from 5.3 eV 
in the atomic limit to 2.6 ± 0.2 eV at a distance of 1.6 

A). 

As we have shown in Sec. Ill B the adsorption of Co 
atoms on graphene differs from the adsorption of closed- 
shell systems because it involves transitions between dif- 
ferent electronic configurations. The nature of these tran- 
sitions can be interpreted in terms of the two primary 
effects: (i) Coulomb repulsion between the electrons in 
the 3c? shell of Co and (ii) electrostatic (Pauli) repulsion 
between the 4s (Co) and it (graphene) electrons. The 
first effect can be described by the Ud-d parameter and 
its evolution with distance, while the second cannot be 
captured by the model Hamiltonian [Eq. (|3J)] since it does 
not explicitly contain the surface energy levels. Since the 
4s-7r orbital repulsion is essentially an electrostatic ef- 
fect, it can be characterized by the 4s-7r orbital splitting 
obtained by means of the simple HF approach Jfi^r^ 
whose variation with distance is also shown in Fig. [6] 

Analyzing the behavior of Ud-d and 5 S - V , one can see 
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that both effects tend to move Co efectrons from the 4s to 
the 3d shell, which leads to the interconfigurational tran- 
sitions of the TM atom upon adsorption. At large dis- 
tances, the variation of the two factors is small, and the 
most favorable configuration (3c? 7 4s 2 ) is governed by the 
strong repulsion in the 3c? shell. As the atom comes closer 
to the surface, the occupation of the 3d shell becomes 
more favorable due to the reduction of the Coulomb re- 
pulsion, whereas the occupation of the 4s shell becomes 
less favorable due to the increasing Pauli repulsion be- 
tween the 4s and 7r electrons. This gives rise to the 3o? 8 4s 1 
and 3c? 9 configurations of Co. 

It is important to emphasize that both described effects 
have a non-negligible contribution and play an important 
role in the determination of the ground-state electronic 
configuration of the adsorbate. In comparison to the elec- 
trostatic repulsion between Co and graphene, the screen- 
ing of the Coulomb repulsion is essentially a correlation 
effect and thus represents a challenge for one-electron 
theories. In the next section, we analyze the perspec- 
tives of the DFT-based approaches in the description of 
TM atom adsorption on graphcnclikc structures. 



D. Connection to DFT 

As a final step, it is interesting to draw a parallel be- 
tween the results obtained in the present work and the re- 
sults of previously reported DFT studies. The standard 
GGA calculations predict that the 3d 9 configuration is 
the ground state of Co on graphene^^ The same con- 
figuration persists if a small [/-correction (U = 2 eV) 
is applied within the GGA+U method^ The applica- 
tion of larger corrections ({7 = 4 eV) yields a configura- 
tion close to 3c? 8 4s 1 .— ^ Finally, hybrid- functional calcu- 
lations, which incorporate some fraction of the Hartree- 
Fock exchange energy, also predict the 3g? 8 4s 1 configura- 
tion for Co on graphene.— 

As one can see, the electronic configurations of Co ob- 
tained within the DFT for small and intermediate values 
of U (0-4 eV) are consistent with the configurations that 
correspond to the same values of U in the present work 
(see Fig. [5] and the discussion above) . This agreement 
justifies the use of on-site Coulomb corrections in DFT 
for the description of TM atoms on graphene. However, 
being significantly dependent on U, the DFT results are 
applicable only to a certain range of adsorption distances, 
while the entire adsorption energy curve cannot be re- 
produced in a consistent way<^ This should be carefully 
taken into account in further studies, especially in the 
determination of ground-state configurations. 

The inability to describe the total potential energy 
curves for TM atoms on graphene within the DFT also 
leads to difficulties in the determination of the binding 
energies. Since an energy calculation of an isolated TM 
atom requires its own U correction to be reliable, the 
comparison of energies [e.g., by using Eq. ([T])] becomes 
meaningless^ On the other hand, if the same U is used 



for the calculation of energies of the interacting and non- 
interacting systems, the binding energies are strongly 
overestimated due to an improper description of isolated 
TM atoms. Another point, though less important, that 
prevents DFT from giving accurate binding energies is 
the lack of nonlocal dispersion interactions in the stan- 
dard DFT approximations (LDA/GGA). 

In summary, the deficiency of the DFT-based ap- 
proaches with respect to the description of TM atoms 
supported on graphenelike structures can be attributed 
to their inability to handle electron correlations in a con- 
sistent and systematic way. 



V. CONCLUSION 

In this work, we have applied multireference quantum 
chemical methods to the systematic investigation of ad- 
sorption of a single cobalt atom on graphene. It is found 
that the electron correlation effects play an important 
role in the adsorption. Far from the surface, these effects 
are not significant and manifest themselves only in the 
weak attractive interactions of the van der Waals type. 
In the vicinity of the surface, the role of correlations is 
more pronounced and can be attributed mainly to the 
screening of the strong Coulomb repulsion between elec- 
trons in the 3g? shell of Co. Due to the significant reduc- 
tion of the Coulomb repulsion, the Co states with higher 
3g? shell occupancy become more favorable at short dis- 
tances. As a result, the total adsorption energy curve of 
Co on graphene is not monotonic and exhibits two min- 
ima, which correspond to the 3g? 7 4s 2 (high-spin) and 3c? 9 
(low-spin) configurations. Therefore, two stable states of 
cobalt on graphene are predicted even for a single (hol- 
low) adsorption site. 

We have shown that the calculated adsorption ener- 
gies for Co/graphene are much smaller (i7 a d s < 0.35 
eV) than was reported previously by DFT studies. The 
main reason is associated with the inability to properly 
describe the total adsorption energy curve for TM ad- 
sorbatcs within the DFT. This, in turn, is related to 
the inconsistent treatment of the strong correlation ef- 
fects, even in the presence of on-site Coulomb corrections 
(DFT+{7 methods). 

Although only the case of the Co adatom is considered 
in this work, we believe that the importance of the corre- 
lation effects is not specific only to the Co adsorption. We 
expect qualitatively similar adsorption behavior for other 
TM atoms on graphene, particularly with respect to the 
interconfigurational transitions and relatively small bind- 
ing energies. The determination of absolute values as well 
as the ground-state electronic configurations of different 
TM adatoms on graphene remains an open issue for fur- 
ther studies. 



11 



VI. ACKNOWLEDGMENTS 

The authors thank Tim Wehling, Vladimir Mazurcnko, 
and Danil Boukhvalov for helpful discussions. Support 
from the DFG Priority Program No. 1459 "Graphene" 
(Germany), from Stichting voor Fundamenteel Onder- 
zoek der Materie (FOM, the Netherlands), and from the 
Russian Scientific Program No. 12.740.11.0026 is grate- 
fully acknowledged. 

Appendix A: Determination of electron-electron 
interaction parameters 

According to the Hamiltonian given by Eq. ©, the 
energies of the five cobalt states of interest can be written 
as 

Eflz/l = 2e Al + 2e El + 3e E2 + 2e 4s 

+2U d - d + UJ d - d + 7J s - d + 19U' d _ d , (Al) 



E S=3/2 = e M + + 3£e 2 + £ 4s 

+W d - d + 13J d - d + 5 J s _ d + 25U d _ d , (A2) 



Eflt/l = 1e Al + 3e El + Ae E2 

+4U d _ d + 16J d _ d + 32U' d _ d , (A3) 



^1=1/2 = ^ £ Ai + 2e El + 3e E2 + e 4s 

+2U d - d + 9J d - d + 7J s - d + l9U' d _ d , (A4) 

^1=1/2 = £a -i + ^ e E x + 3e_B 2 + Sis 

+3U d _ d + 13 J d _ d + 3J s _ d + 2W' d _ d . (A5) 

Substracting Eq. ([AT]) from Eq. (|A4]) . and Eq.(|A2]) 
from Eq. (|A5[) . an expression for the exchange param- 
eters can be readily obtained: 

J d - d = \(Ef:tf 2 -Efltf 2 ), (A6) 

Js-d = 2 ( E S d =3/2 ~ E S=l/2) ■ ( A7 ) 

Summing up Eqs. ([AT]) and (|X3)l . subtracting Eq. (|X2|) 
twice, and exploiting the rotation invariance condition 
(U d _ d = U d - d — 2J d - d ), one obtain an expression for the 
intrasite Coulomb repulsion parameter: 

7-r _ rpZifAs 2 , j7T3d 9 4s° 7T n 3ci 8 4s 1 
Ud-d — ^S=3/2 + ^S=l/2 _ ^-^5=3/2 

-2A.A t - El - A E2 - El + J d - d + iJ s - d , (A8) 
where Ai_j = Ei — Sj are crystal-field splitting parame- 
ters. Therefore, assuming that the rotation invariance 
condition is valid for the system under investigation, 
the electron-electron interaction parameters can be ex- 
tracted from the total energies of different states up to 
the crystal-field splitting parameters. 
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